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Abstract 

A novel Quantum Similarity Measure (QSM) is constructed, based on concepts from Information 
Theory. In an application of QSM to atoms, the new QSM and its corresponding Quantum 
Similarity Index (QSI) are evaluated throughout the periodic table, using the atomic electron 
densities and shape functions calculated in the Hartree-Fock approximation. The periodicity of 
Mendeleev's Table is regained for the first time through the evaluation of a QSM. Evaluation of 
the Information Theory based QSI demonstrates however that the patterns of periodicity are lost 
due to the renormalization of the QSM, yielding chemically less appealing results for the QSI. A 
comparison of the information content of a given atom on top of a group with the information 
content of the elements in the subsequent rows reveals another periodicity pattern. 

Relativistic effects on the electronic density functions of atoms are investigated. Their impor- 
tance is quantified in a QSI study by comparing for each atom, the density functions evaluated 
in the Hartree-Fock and Dirac-Fock approximations. The smooth decreasing of the relevant QSI 
along the periodic table illustrates in a quantitative way the increase of relativistic corrections with 
the nuclear charge. 

PACS numbers: 
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I. INTRODUCTION 



Quantification of molecular similarity using an electron density based Quantum Similar- 
ity Index (QSI) yields compact information on the similarity in shape and extent of the 
electron density distribution of various molecules. These data can be used as descriptor in 
comparative discussions of molecular structure and reactivity [l| . In recent years one notices 
a multitude of methodological studies on various aspects of Quantum Similarity of molecules 
such as the use of different separation operators {2}, the replacement of the density by more 
appropriate, reactivity oriented functions jjj, 0] , within the context of conceptual DFT [s| 
and the peculiarities arising in the comparison of enantiomers js, 3]. 

Most remarkably the field of similarity of isolated atoms remained practically unexplored, 
with the exception of two papers by Carbo and coworkers jfcj, ^ and a third one by the present 
authors [10j . In the first paper Q| atomic self-similarity was studied, whereas the second 
one [9| contains a relatively short study on atomic and nuclear similarities, leading to the 
conclusion that atoms bear the highest resemblance to their neighbors in the Periodic Table. 



In 



10| we reported results of the Carbo index, which turns out to mask the periodicity of 



Mendeleev's Table, followed by results of an Information Theory based approach, where 
the Information Discrimination was evaluated for atomic electron density functions. Hereby 
the periodicity in the atomic electron density and shape functions throughout Mendeleev's 
Table was regained. The relevance of Information Theory in applications to electron density 
functions of atoms has also been confirmed in a study of the iV-derivative of the Shannon 

n 

entropy of shape functions [11] . 

The present work is, in a first part, concerned with the construction of a new Quantum 
Similarity Measure (QSM) along the lines of the study in 10j. The simplification of the 
density function based Quantum Similarity Index to a shape based expression emphasizes 
the potential of the shape function as an alternative to the density function. The use of 
the shape function as a fundamental descriptor of atomic and molecular systems is indeed a 
current topic of investigation in the domain of Quantum Chemical description of atoms and 

nnfl 

molecules [12l . 1 131 Il4j |. After defining the mentioned Information Theory based QSM and 
the corresponding QSI, they are evaluated for all pairs of atoms in the Periodic Table. The 
results are interpreted and investigated for patterns of periodicity. 

In the second part of this work we investigate the relativistic effects on the electron 
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density functions of atoms and their quantification using QSI. From the relativistic effects 
on total ener gies one can infer these effects have implications for the electron densities, as 
visualized in llj for the Rn atom. The effect of relativity on atomic wave functions has 



been studied in the pioneering work of Burke and Grant [15| who presented graphs and 
tables to show the order of magnitude of corrections to the hydrogenic charge distributions 
for Z = 80. The relative changes in the binding energies and expectation values of r due to 
relativistic effects are known from the comparison of the results obtained by solving both 
the Schrdinger and Dirac equations for the same Coulomb potential. The contraction of 
the ns-orbitals is a well known example of these relativistic effects. But as pointed out by 
Desclaux in his u Tour historique" [161 ]. for a many-electron system, the self-consistent field 
effects change this simple picture quite significantly. Indeed, contrary to the single electron 
solution of the Dirac equation showing mainly the mass variation with velocity, a Dirac- Fock 
calculation includes the changes in the spatial charge distribution of the electrons induced by 
the self-consistent field. The framework of QSI offers a simple way of quantifying relativistic 
effects on atomic electron densities via a comparison between non relativistic Hartree-Fock 
and relativistic Dirac-Fock electron density functions. 

II. METHODOLOGY 

A. Quantum Similarity Indices and Information Theory 

Our work is situated in the context of a mathematically rigorous theory of Quantum 
Similarity Measures (QSM) and Quantum Similarity Indices (QSI) as developed by Carbo 
{2I, [n]]. This theory encompasses quantum objects, e.g. atomic and molecular systems. 
Following Carbo, we define the similarity of two atoms (a and b) as a QSM Z ab (Q), 

Z ab (Vt) = J p (ri) J](ri, r 2 ) p 6 (r 2 ) dridr 2 (1) 

where f2(ri,r 2 ) is a positive definite operator. Renormalization to 

SIq = /z Wz rn? (2) 
yields a QSI Sin with values comprised between and 1. 

By choosing the operator f2(rx,r 2 ) to be the Dirac- 5(r 1 -r 2 ) function, expression [1] reduces 
to an overlap integral, yielding the simplest form of the Carbo Similarity Index, (Sis), 



after normalization. Using shape functions defined as a = p/N (N being the number of 
electrons), the QSI for density functions simplifies to a QSI for shape functions. This yields 
the important result that by investigating the similarity between two systems, we are in fact 
comparing their shape functions. This motivates the investigation of the shape based QSM , 
found by substitution of p by a in expression [TJ 

For the construction of a new QSM, we considered the introduction of concepts from In- 



formation Theory 
research in general 



19 L which has recently been of increasing relevance to quantum chemical 
and to the investigation of the electron densities in position and 



momentum space in particular. (A thorough discussion can be found in the pioneering work 
of Gadre and Sears [20(.) In our previous work we reported on the information entropy of 
atomic density and shape functions respectively defined as 

AS>= [ Pair) log J^-dr (3) 



% = I a a {v)\og^\dv (4) 
J oo W 



with po( r ) the density of the prior or reference atom. As motivated in 10] we set the 
density function of the prior equal to the density of the noble gas preceding the atom under 
investigation in the periodic table, scaled by the factor where N a and N are the number 
of electrons, respectively of atom a and its reference. In this way the prior density po( r ) and 
the density of atom a, p a (r), yield the same number of electrons upon integration. It was 



shown [10| that these quantities reflect the periodic evolution of chemical properties in the 
Periodic Table and that Kullback's interpretation can be formulated in terms of chemical 
information stored in the density functions when we make this particular choice for the prior 
densities. 

n 

Following the conclusions in [lQ|, one can see that it would be interesting to compare the 
Information Entropy, evaluated locally, AS^(r) = p a (r) log , of two atoms by use of a 

QSM. To that purpose the integrand in expression [3] is considered as a function, which gives 
the Information Entropy locally at a given point r. The construction of the corresponding 
QSM becomes straightforward by considering the overlap integral (with Dirac 5 as separation 
operator) of the local Information Entropies of two atoms a and b 

Z ah {5) = [ p a (r)log 1 ^^p fe (r)log 1 |^Trfr= f AS£(r)ASf(r)dr. (5) 



A QSI can be denned by normalizing the QSM as before, via expression [2j The QSM and 
the normalized QSI give a quantitative way of studying the resemblance in the information 
carried by the valence electrons of two atoms. Expression can be rewritten in the form [2] 
by identification of the operator fi[p a (ri), p&(r 2 ); p (ri), Po( r 2)] = In Z a{r f > In - 5(r a - 

r 2 ), where we explicitly write the functional dependence on p a ( r i) and p&(r 2 ) and the para- 
metrical dependence on po( r i) and Po(r 2 ). 

The obtained QSI trivially simplifies to a shape based expression 

/ASJ(r)AS?(r)dr 



5/, 



(5) 



(6) 



7 AS?(r) ASz(r)dTyJf A^(r) A^(r)dr 

where a shorthand notation is used by omitting the explicit dependency of r. The simplifica- 
tion can be generalized from the local information distance operator AS p (r) to any operator 
F p (r), which is linear in p (thus satisfying F p (r) = NF a (r)), as follows : 

/ FZ(r)F b p (r)dx j F°{r)F h °{v)dr 



J F p (r) Ft{r)dr J J F b p (v) F p {v)dv J J F*(r) F° a {v)dr J J F b °(v) F b °(v)dv 



(7) 



In agreement with the fact that the shape function completely determines the properties 
of a system, as discussed in [2l|, the relevance of the QSI as a tool to compare physical 
properties of atomic electron density functions is confirmed. This characteristic distinguishes 
the QSI above, together with the Carbo QSI from other similarity measures (e.g. Euclidian 
distance [22j|, Tanimoto [231] and Ho dgkin- Richards [241]). 



B. Atomic electron density functions 

1. Non-relativistic atomic electron density functions 

The atomic electron density functions were evaluated from non-relativistic numerical 
Hartree-Fock wave functions optimized on the L^-term ground state of neutral atoms (nu- 
clear charge 3 < Z < 103), as specified in the Table of electron configuration and term 



value given by Bransden and Joachain or the NIST website [25|, |26[. The extension of 
the original Froese-Fischer's code by Gaigalas [27| allows the calculation of term-dependent 
Hartree-Fock orbitals for any single open subshell case. However, for the two-open subshells 
cases, this version is limited to (ns)(n'l) N (I = 0, 1, 2, 3, . . .), (np) N (n r l) (I = 0, 1, 2, 3, . . .) 
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and (nf)(n'd) configurations. This computer code then covers all the ground levels of the 
periodic table, except atoms with ground configuration f N - 2 d. For these specific systems, 
we used the "MCHF atomic-structure package ATSP2K" 281 ]. relying on the combination 
of the second-quantization approach in the coupled tensorial form, the generalized graphical 
technique and angular momentum theory in orbital, spin and quasispin spaces and on the 
use of reduced coefficients of fractional parentage [291 ] . 

In the non-relativistic Hartree-Fock approximation, in its single-configuration version, the 
atomic wave function is limited to one configuration state function (CSF) [30I 

\anLSM L M s ), 

simultaneous eigenfunction of the inversion operator, the orbital angular momentum L 2 , the 
spin angular momentum S 2 and their projections L z and S z , that can be built from the 
one-electron spin-orbitals 

V<(r) = R nl {r)Yr(6,ct>)Xi/2,m a (<7) = rP„^>T(M)Xi/2, m >). (8) 

a denotes all the information needed to specify unambiguously the term considered (config- 
uration and coupling tree). The optimized one-electron numerical radial functions {P n i{f)} 
are used to determine the corresponding LS'-dependent electron density function from the 
following expression 

»M = I £ ^ (9) 

nl 

where q n i is the occupation number of the subshell considered. In the case of uncompletely 
filled subshells, spherical averaging over the (Ml, Ms) term components was applied, yielding 
a spherical electron density function, as elaborated in jlO[ . In an LS'-dependent Hartree Fock 
scheme, the radial wave functions are allowed to vary, for a given electronic configuration, 
from one term to another. Eq. [9] should then strictly be read as 

p(r) = PaLS (r) (10) 

An easy way of testing the calculated density function p(r) is to check that its integration 
yields the total number of electrons 

POO 

4tt / p(r)r 2 dr = J2 ( lm = N, (11) 

J ° nl 

as expected from the normalization constraint of the Hartree-Fock numerical one-electron 
radial wave functions. 
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2. Relativistic atomic electron density functions 



For the purpose of quantifying the relativistic effects on the electron density functions, 
we evaluate the similarity of Hartree-Fock and Dirac-Fock density functions using a point 
nucleus aproximation. In the relativistic scheme, the atomic wave function is, in the most 
general combination of configuration state functions (CSF's) 

\nJMj) = ^c v \utiJMj) (12) 

eigenfunction of the inversion operator, the total angular momentum J 2 and its projection 
J z . v denotes all the necessary information for specifying the relativistic configuration 
completely. The CSF are built on the one-electron Dirac four-spinor 

, , x 1 ( Wx£(n) , 

^( r ) = - , , „ * ( 13 ) 



r 



where is a two-dimensional vector harmonic. It has the property that Kipi(r) = 

Kipi(r) where K = f3(cr ■ L + 1). 

The large {P(r)} and small {Q{r)} components are solutions of a set of coupled integro- 



differential equations 3jJ. The mixing coefficients {c u } are obtained by diagonalizing the 



matrix of the no-pair Hamiltonian containing the magnetic and retardation terms [32j . The 
two coupled variational problems are solved iteratively. For a complete discussion on rela- 
tivistic atomic structure we refer to (3^]. The present calculations have been performed using 



the MDF/GME program of Desclaux and Indelicato [34j including both the magnetic and 
retardation part of the Breit interaction in the self-consistent process, but not the vacuum 
polarization. 

It is to be noted that the relativistic scheme rapidly becomes more complicated than 
the corresponding non-relativistic one. For example, if the ground term of Carbon atom is 
described, in the non-relativistic one-configuration Hartree-Fock approximation, by a sin- 
gle CSF \ls 2 2s 2 2p 2 3 P), the relativistic equivalent implies the specification of the J- value. 
For J = corresponding to the ground level of Carbon, the following two-configuration 
description becomes necessary 

\ u ls 2 2s 2 2p 2 " (J = 0)) = c 1 \ls 2 2s 2 (2p*) 2 (J = 0)) + c 2 \ls 2 2s 2 2p 2 (J = 0)), 
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implicitly taking into account the relativistic mixing of the two LS-terms ( 1 S and 3 P) 
arising from the 2p 2 configuration and belonging to the J = subspace. p* and p correspond 
to the j-values, j = 1/2 (k = +1) and j = 3/2 (re = —2), respectively. 
By averaging the sublevel densities, 

'M = pTTT) £ » JMJ ^ < 14) 

V > Mj=-J 

the total electron density becomes spherical for any open-shell system, as found in the non- 
relativistic scheme lCj, and can be calculated from 



^ = ^3^^ (15) 

where q nK is the occupation number of the relativistic subshell (nn). 

III. INFORMATION THEORY QSM AND QSI FOR ATOMIC DENSITY FUNC- 
TIONS 

In this section the results of the QSM and QSI, evaluated for shape functions of all pairs 
of atoms in the periodic table are discussed. To facilitate the interpretation of the results 
of the Information Theory based QSM and QSI a graphical representation of the Carbo 
QSM (figured]) and QSI (figure [2]), already mentioned in [h]], is given. The results of all 
the possible pairs of atomic shape functions are given in a 3 dimensional graph, where the 
vertical axis indicates the QSM or QSI value of the atoms with nuclear charges Z a and Zj,. 
These figures show that the general trend for any fixed Z a is similar for all atoms: all the 
cross-sections of the 3 dimensional graphs show the same evolution. In figure [3] we show 
the cross-section of the results of the QSM between Pb [Z = 82) and all other atoms. The 
general trend of the overlap QSM of the density functions increases monotonically with 
increasing volume of the atoms, as pointed out in {9]. 

The evaluation of the Information Theory based QSM (figure H]) was found to be positive 
for all investigated atoms reveals a picture corresponding to the periodicity of Mendeleev's 
Table, which can be distinguished by looking at the cross-section for Pb in figure The 
results correspond to the evolution of chemical properties first of all in the sense that for 



each period the QSM increases gradually from the first column to the last, 
and Hardness are properties which reveal a similar evolution throughout 
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FIG. 1: Overlap integral of the non-relativistic Hartree-Fock shape functions (the QSM appearing 
in the definition of the Carbo QSI) evaluated for all pairs of atoms. A monotonic trend of increasing 
QSM for heavier atoms is revealed. The vertical axis corresponds to the QSM of the atoms with 
nuclear charges Z a and Z\> given by the axes in the plane. 

the sense that neighboring atoms with large nuclear charge differ less than neighboring light 
atoms, e.g. the difference between the QSM values of two atoms in the first period is large 
in comparison to the difference in QSM between two neighboring Lanthanides. Considering 
all the cross-sections of Figure H] reveals that the periodicity is regained throughout by the 
choice of the reference atoms, as it yields low QSM values for atoms similar to the chosen 
prior. One notes however that the QSM does not reveal results, which reach maxima when 
a given atom is compared with another atom of the same group. Moving to the QSI, the 
periodicity of the QSM is lost due to the normalization (figures |6] and Cj) . In figure [7] the 
change of prior is still visible due to the gaps (discontinuities) at the positions where the 
prior changes, but the normalization blends out the clear periodic evolution of the QSM 
in graph [51 This leads to the conclusion that the normalization, which yielded the nearest 
neighbor effect for the Carbo QSI in figure EJ can overwhelm the characteristics of a QSM. 
Changing the point of view, we can opt to investigate which atom of a given period of 
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QSI 




FIG. 2: The Carbo QSI evaluated for all pairs of atomic non-relativistic Hartree-Fock shape func- 
tions, revealing a nearest neighbour effect The vertical axis corresponds to the QSI of the 
atoms with nuclear charges Z a and Zf,, indicated by the axes in the plane. 



the table belongs to a certain column and in which way the atoms should be ordered within 
the period. This can be done by investigating the QSI with the top atoms of each column 
as prior. Formulated in terms of Kullback Liebler information discrimination the following 
is evaluated. For instance, when we want to investigate the distance of the atoms Al, Si, 
S and CI from the iV-column (group Va), we consider the information theory based QSI 
in expression [51 where the reference densities po and po> are set to pn, Pa to pau Psu Pp, 
etc. respectively and p# to pp, i.e. we compare the information contained in the shape 
function of iV to determine that of P, with its information on the shape function of Al, 
Si, S, CI. The data in table [I] reveal a 1. for the element P (by construction) with values 
continuously decreasing from unity for the elements to the left and to the right of the 
iV-column. This pattern is followed for the periods 3 up to 6, taking As, Sb and Bi as 
reference, with decreasing difference along a given period (see Figure [8]) in accordance with 
the results above. Note that the difference from 1. remains small, due to the effect of the 
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FIG. 3: The cross section of the graph in figure Q] for Zb = 82. The vertical axis gives the QSM 
value for Pb with the atom with nuclear charge Z a . 



renormalization used to obtain the QSI. 

TABLE I: Numerical results of the QSI with prior atoms set to the elements on top of the columns. 
The information present in the shape function of N to obtain information on that of A is compared 
with the information present in the shape function of N to obtain information about iV-group atom 
of the corresponding period. 



Al: 


0.98656 


Si: 


0.99688 


P: 


1. 


S: 


0.99735 


CI: 


0.99031 


Ga: 


0.99880 


Ge: 


0.99971 


As: 


1. 


Se: 


0.99973 


Br: 


0.99897 


In: 


0.99957 


Sn: 


0.99989 


Sb: 


1. 


Te: 


0.99990 


I: 0.99961 


Tl: 


0.99986 


Pb: 


0.99996 


Bi: 


1. 


Po: 


0.99996 


At: 


0.99987 
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FIG. 4: Information Entropy based QSM for pairs of atoms in the periodic table, with the noble gas 
of the previous row as prior for each given atom. A clear periodic character can be distinguished. 
A non-relativistic Hartree-Fock approach was used. 



IV. INVESTIGATION OF RELATIVISTIC EFFECTS 



In this section we discuss the relativistic effects on atomic electron density functions. We 



first illustrate the difference of the radial density functions D(r) defined as 25|, 



D(r) =A7rr 2 p(r), (16) 

calculated in the Hartree-Fock (HF) and Dirac-Fock (DF) approximations for the ground 
state 6p 2 3 P of Pb I (Z — 82) according to equations [9] and [151 respectively. These are 
plotted in figure EE which shows the global relativistic contraction of the shell structure. 
Another interesting observation can be done from figure \TU\ displaying, as was done in 



ll| for the Rn atom, the accumulated difference between the DF and the HF radial density 



functions [T6l denned as 



AD(r) = I ( D DF (r') - D RF \r') ) dr' (17) 



o 
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FIG. 5: The cross section of the graph in figure H] for = 82. The periodic character is regained. 

as a function of r, the radial distance to the nucleus. We see in this way that there is an 
excess charge, varying between and 0.9 due to relativistic effects. One notices that the 
contraction of the total radial density function reveals a shell structure. Since the densities 
are normalized to the same number of electrons, the accumulated difference converges to 
for large values of r. 

Employing the framework of QSI to compare non-relativistic Hartree-Fock electron den- 
sity functions p HF (r) with relativistic Dirac-Fock electron density functions p DF (r) for a 
given atom, the influence of relativistic effects on the total density functions of atoms can 
be quantified via the QSI defined below 



Z HFtDF (8)=fp HF (r)p DF (r)dr (18) 
Sis = , Zhf - df ^ , (19) 

y Zhf,hf{&)\/ Zdf,df(5) 

where 5 is the Dirac-5 operator. 

In figure [TT1 we supply the QSI between atomic densities obtained from numerical Hartree- 
Fock calculation and those obtained from numerical Dirac-Fock calculations, for all atoms 
of the periodic table. The results show practically no relativistic effects on the electron 
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FIG. 6: Information Entropy based QSI for pairs of atoms in the periodic table, with the noble 
gas of the previous row as prior for each given atom. The vertical axis corresponds to the QSI of 
the atoms with nuclear charges Z a and given by the axes in the plane. 

densities for the first periods, the influence becoming comparatively large for heavy atoms. 
To illustrate the evolution through the table the numerical results of the carbon group 
elements are given in table [III and highlighted in the graph in figure [TTJ From the graph it 
is also noticeable that the relativistic effects rapidly gain importance for atoms heavier than 
Pb. 

Investigation of the convergence of the QSI in function of the radius r can shed some 
light on the importance of core region of the density functions for the QSI. In figure [12] we 
plot, for the Pb atom, the numerical results of the QSI defined as 

ZhfM* r) = lo L P HF (*') P DF ( r ') dndr ' ( 20 ) 
SI s (r) = (21) 

y Z HF>HF {&;r)^/ Z DF}DF (S;r) 

where the integration over Q represents the integration over all angles. The plot shows a 
very fast convergence, the total QSI value being reached already for r = 0.2 a . This result 
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FIG. 7: A cross section for Z\> = 82 of the graph of Information entropy QSI in figure [H The 
change of reference atom is still visible, but the periodicity is not regained. 



TABLE II: Numerical results of the QSI for the carbon group, highlighted in graph [TTJ 



Atom 


QSI 


C 


0.99999 


Si 


0.99996 


Ge 


0.99945 


Sn 


0.99661 


Pb 


0.96776 



demonstrates the dominance of the inner region of the density function for this type of QSI. 
This picture does not reflect the influence of relativistic effects on the valence electrons, 
which is visible in the accumulated difference in figure El The rapid convergence of the QSI 
can be accounted to the fact that the densities, in the overlap integral, themselves are much 
larger in regions of small radius r, whereas the clear influence of relativistic effects in the 
accumulated difference picture is due to the fact that the difference between Hartree-Fock 
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FIG. 8: Results of the information theory based QSI with the atom on top of the column as prior. 
The symbol in the legend indicates the period of the investigated atom and the nuclear charge 
Z-axis indicates the column of the investigated atom. (For example Ga can be found as a square 
Z = 5). 

and Dirac-Fock radial densities remains of the same order of magnitude, converging to zero, 
as can be seen in figure [9j 

V. CONCLUSION 

In this communication we report on the development and calculation of a new, Informa- 
tion Theory based, Quantum Similarity Measure (QSM) and its corresponding Quantum 
Similarity Index (QSI) for atoms, using their density functions and shape functions. We 
show that a QSM constructed with the Kullback Leibler Information Entropy loses its peri- 
odic character upon normalization. One might say that the normalisation renders the QSI 
insensitive to certain characteristics, clearly present in the QSM. To regain the periodicity 
with the information theory based QSM, the choice of the prior for each atom as the density 
of the noble gas of the previous row, normalized to the same number of electrons in the 
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FIG. 9: DF and HF density distributions D(r) = 47rr 2 p(r) for the neutral Pb atom (Z = 82). The 
contraction of the first shells is clearly visible. 
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FIG. 10: Accumulated difference between the DF and the HF densities of Pb. 
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FIG. 11: Similarity of non-relativistic Hartree-Fock with relativistic Dirac-Fock atomic density 
functions with highlighted results for the CIV group atoms. 



atom under investigation, is crucial. The results of the QSM are in agreement with chemical 
intuition in the sense that the difference in QSM of two successive light atoms is large in 
comparison to the difference in QSM of two successive heavy atoms, meaning that light 
atoms next to each other in the Table differ more than neighboring heavy atoms. When 
looking at the results of Lanthanides and Actinides in particular we find high similarities 
indeed. This interpretation is not regained by looking at the QSI, with the prior set to the 
noble gas atoms. It is rewarding that the comparison of information content of the shape 
function of a given top atom in a column with the atoms of the subsequent period(s) reveals 
another periodicity pattern. 

The visualization of the influence of relativistic effects on the radial density distribution 
reveals a shell structured excess charge, corresponding to the contraction of the charge 
distribution. The importance of relativistic effects for the electron density functions and 
shape functions of atoms has been quantified via a study based on QSI. A plot of the QSI as 
a function of the nuclear charge shows that the densities of light atoms are barely influenced 
by including the relativistic corrections and that the influence of relativistic effects increases 
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FIG. 12: Convergence of the QSI of HF and DF densities of Pb, as a function of r, revealing the 
large influence of the inner region of the density function. 

monotonically with the nuclear charge of the neutral atom throughout the Periodic Table. 

As suggested by the investigation of: i) the Carbo QSI, ii) the Information Theory based 
QSI and iii) relativistic effects via the Carbo QSI for atoms, the Carbo QSI reflects the 
similarity of the core region of the density function, i.e. it fails to reflect the importance 
of the valence electrons, which is essential from a chemical point of view. It would be 
interesting to investigate if the valence region can be given more weight in a similarity 
study, by introducing an appropriate separation operator in the definition of the QSI. 

Although correlation effects are neglected in the present work, a similar QSM/QSI ap- 
proach can be used for investigating how much the electron densities are affected by cor- 
relation, comparing the same atom in the single- and multi-configuration non relativistic 
Hartree Fock approximations. 
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